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Abstract 

A varying external field approach in QCD sum rules is formulated in a sys¬ 
tematic way to treat the weak decay form factors and their dependence 
in the process of i? —> p t From the form factor sum rules, we can also 
obtain the mass sum rules for B and p mesons, which can help us determine 
the reliable Borel windows in studying the relevant form factor sum rules. 
In this way, we thus demonstrate that some QCD sum rule calculations in 
the literature are less reliable. We also include induced condensate contri¬ 
butions, which have been ignored, into the relevant sum rules. We obtain 
the ratios r(.B° ^ p'^e~i>e)/ ^ 7r+e“r'e)~ 0.94 and T{B^ p'^T~Vr)/ 

~ 1.15. We apply this approach to re-examine the case of 
the D meson decay. 
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I. INTRODUCTION 


Semileptonic decays, because of their simplicity, provide an excellent laboratory where 
physicists can study the effect of nonperturbative QCD interactions on the weak decay 
process. A detailed understanding of these processes is also essential for determining the 
magnitudes of CKM quark-mixing matrix elements. 

The weak decay form factors of B ^ p i have been calculated by using various 
approaches. However, the results obtained from the traditional sum rules |]^ by considering a 
three-point correlation function with suitable interpolating current seem to conflict seriously 
with others’ theoretic results, such as light-cone sum rules [|],|1, lattice simulations @], quark 
models [^-P], and the external field approach of QCD sum rules |^. Recently, Ball and 


Braun [|^] re-examined such a process by studying the light-cone sum rules with the modihed 
p meson wave functions. In their study, all of the soft (non-perturbative) parts are absorbed 
into the p meson wave functions. The accuracy of their results is dependent on the shape of 
the wave functions. Similarly, within the pQCD approach |^, the light meson is described 
by a phenomenological model function which can be taken, e.g., from QCD sum rules. 

In this work, we shall use the varying external held approach of QCD sum rules which 


has been developed earlier []^ to re-examine the weak decay form factors for B —>■ p i ui 
with a complete calculation. This approach, in spirit, is similar to Ref. [|^. We will present 
the idea of the induced condensates and show how the induced condensates enter the sum 
rules. In the calculation we hnd an additional contribution from the induced condensates, 
which was ignored in Ref. |^. Consequently, the Ai form factor will obtain an additional 
10% contribution from the induced condensate, while the other form factors will not. In 
addition, we also propose a reliable method that can be used to determine the Borel windows 
in studying the relevant form factor sum rules. 

The rest of this paper is organized as follows. The concept of the induced nonlocal 
condensates will be formulated in the part (a) of Sec. II. In the part (b) of Sec. Ill, the 
external held approach of QCD sum rules will be set up to investigate the various decay 
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form factors of B ^ p i ug. In Sec. IV numerical results and discussion are presented. Sec. V 
contains a brief summary. 


II. THE METHOD 


A. Quark propagator in the presence of an external variable field 


The quark propagator in the external vector (axial vector) held is described by the 
additional term, in the Lagrangian. 


AC\x) = 


( 1 ) 


or 


AC\x) = 


( 2 ) 


where 


Ja = U^ab, Ja = U'^albb- 


( 3 ) 


Here V" {A^) and q°‘ are the amplitude and momentum of the external vector (axial-vector) 
held, respectively, b stands for the b quark held operator, and u is for the u quark held 
operator. Hence the quark propagator depicted in Fig. 1(a) can be written down as 


(i)Syo)>y,^, = (T<(i)65(0))''% + {:(i)Syo) (4) 


with the notation (...) = (0|... |0). Here a and b are the color indices while a and (3 the 
Lorentz indices. The hrst term depicted in Fig. 1(b), on the right hand side of Eq. (4), 
can be calculated from perturbation theory, while the second term depicted in Fig. 1(c) is 
the induced condensate dehned through this paper. Neglecting the radiative corrections, we 
obtain 


(T<(a;)6;(0))r = 


- ml)[{p + qY - ml] 


( 5 ) 
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Using the identity 


= l(bu) + 75 ( 675 ^) + Yihpu) 
= l{bu) + J5(b'y5u) + Yihpu) 


Ylb{blplbu) + -aP'"{bapru) 


!ab 


- 0/3 


- Yl^ib-ip-i^u) + -(T^'^b{h(Tpr-i^u) 


zab 


- 0/3 


( 6 ) 


we now project out Eq. (5) on the Dirac matrix basis 1, 75 , 7 p, 7 p 75 , cTp^. Thus the induced 
condensate (the second term on the right hand side of Eq. (^), with the aid of Eq. (^, is 

— — 

(: ni(x)hl(0) O’" = - ^ ((^ KO)«(l) :)i„^ + (: 5(0)75«(j^) 

+ (: K^)lpu{x) :)(7^)o/3 - (: &(0)7p75M(a;) 

1 \ induced 

+ ^{■H^)(^pMx) :){aP^)^pj^ , (7) 

or 

- - - 
{: = - 15 - ((^ K0)75f (x) :>(75).,J 

+ (: H^hpuix) ■){Y)a^ - {■ K^hp'yMY ■){Yi^)ap 

Y \ induced 

+ -(: 6(0)crp^75M(a:) :)(cr^^75)o/3) , (8) 

2 /a 

where 


{: 6(0)r«(x) :>“77 = {T6(0)r«(x))v(7) - (T6(0)r«(x))7j, (9) 


with 


(T6(0)rM(a:))v = J Yzd'i^T{u{zhpb{z),b{0)Tu{x))), 

{TmTu{x))A = -lA^ j dSe*''^(T(h(z)7p75K^),K0)rM(a;))), (10) 


and 


{Tb{0)Tu{x)yY^ = -iN, 


Yp -ipx Tr [(i6 + m„) p{i)+ ^ + mj,)r] 
(27r)4^ {p^ - ml)[{p + qf - ml] 


( 11 ) 


(T6(0)rM(x))^^’'* = -%N, 


Yp ^^p^ Tr[(j6 + m^) A75((i^+ /l) + mk)V] 

(2nY (p 2 - ml)[ip + g )2 - ml] 


( 12 ) 
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Here T is the generic notation for the Dirac matrix basis 1, 75, 7^, 7^75, cxpr and the path- 
ordered gauge factor is implied by 

Q{0,x) = Pexp[—igs [ dax'^AJax)]. (13) 

Jo 

In the case of the hxed-point gauge (the Fock-Schwinger gauge) x^A^[x) = 0, this factor is 
equal to unity. 


B. The Derivation of QCD Sum Rules 


For the form factors of B ^ p i we consider the following two point Green’s function 
in the external variable vector held V or axial-vector held A: 


n;;7(p,p'; = Ad\ j=(0)})v,j), 


(14) 


where (...) = (0|... |0), = d'y^u and = hy^'y^d. The interaction with the external 

vector held is described by the additional term, A£, in the Lagrangian as shown in Eqs. (1) 
and (2). In the following calculation we consider only the amplitude of linear in the 

external held V(^). Generally speaking, the coefficients of tensor structures H* in can 
be expressed as a double dispersion relation form: 

IfApip']q^) = [ ds [ ds'-. — ^ ) _|_ s^htraction terms, (15) 

J J [s — p^){s'— p’^j 

where the subtraction terms have the form 

'A{s')ds' , 2 \ fA{s)ds 


subtraction terms = Pi^Up^) J ^ + P'2i,v(p'‘) J A ^ ^ -Psaclpt?”)- 


(16) 


P'^ J s — p^ 

One hnds that such a contribution from the subtraction terms is of no importance since, 
after performing the double Borel transform, all of the subtraction terms, Pi^u, P 211 U, and 
P^^u) disappear. (For further discussions, see Ref. [|^]) The double Borel transform Rpj] 
on Hj in both variables p^ and gives 


R[n,]= r ds r ds' p, 

Jo Jo 


( 17 ) 
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The hadronic side of Eq. (14) is represented as a sum over the hadronic states. If the Borel 
masses are chosen properly, the hadronic side of Eq. (14) will be dominated by the lowest 
hadronic states with the contributions from the higher states and the continuum suppressed. 
At the hadron level (the right hand side), the spectral function of can be written as 

jBTn? 

- -[f9a^lPu + + p')aPu + a.p^QaPu 

9p 

+other terms ...] x 6{s — m%)5{s' — rn?p) + higher states, 

= V^l^I!^[2ige^pppp^p'f^pP 
9p 

+other terms ...] x 5(s — m%)5{s' — m^p) + higher states, (18) 

where the higher states start from a higher enough value, Sq or Sq. The contribution from 
the higher states could be approximated by the perturbative part, (which starts from the 
value, Sq or Sg) in the OPE of the quark level. The “other terms” in Eq. (18) are the 
irrelevant tensor structures that are not taken into account here and thus are suppressed. 
Here we have adopted the dehnitions 

9 

Tft 

{^\jp\p{pA)) = e(A)M, 

9p 

{p{p, \)\J^ - Ja\B{p)) = 2igeo,pppp'^^p^e\{^ 

-fe*x)a - a+{e*x) ■ P){P + P')a - a.{e*^y p)qa, (19) 

with 

/ = (mp + mB)Ai, 

mp + niB 

a_ = 2mpA, 

2^ _ ^o(g^) - 4l3(g^) . 2^ _ + A2{q^){mB - TUp) 

^ ’ 3(? ) 2mp 

V 

3 = -^- 

"ip + "is 

in the BSW parametrization |Q. 

At the quark-gluon level, the Eq. (14) can be alternatively written as 
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n;;<""(p.p';<(") = -"/e‘'''"{Tr{Sf(0,l)7„SS'''"‘>(a:, 0 ) 7 . 75 }), 


( 20 ) 


where 10,15 


SS(a)(0.^) = / 


d*p , 


(27r) 


ipx 




md 4 2 


(p2 - m2 )2 


- A - 9sX -- 

4 p-irid 


r 


- rud 


V 


+ : h“(0)h^"(0) : +a;^ : <(0)(Zl'^(i^"(0)) : 

+ : ci“(0)(Zl'^Zl^Zl^(^(0)) : 

3! 




2 ! 


( 21 ) 


and 0) = (TM^(a;)6“(0))v(yt) is shown in Eq. (4) but also in the background gluon 

held. Here we have used the hxed-point gauge, = 0 for the gluon held and 

_D" = 


The Feynman integral for the bare loop can be written as 




n;i(p,p';g^) 


3iV“ 

3U“ 


d^k Tr[(/+ /; + mu)ja{^+ /: + mb)7^75(/; + rudh^] 
(27r)'^ k‘^{k + p'y[{k + pY — ml] 

d^k Tr[(/+ /c + m^) 7 ct 75 (i>+ /c + mb)'yu75iA + ^d)7/.] 

(27r)"^ _|_ j 9 /) 2 J('/. _|_ PP _ .^ 2 j 


( 22 ) 


where and have been neglected. The relevant diagram is shown in Fig. 2(a). Using 
the Cutkosky’s rule |^, the integral can be solved easily by taking the imaginary part of 
the quark propagators: l/(p^ — ml + ie) —>■ —2'Kid{p^ — ml). Thus this integral may be 
recast as a double dispersion relation 


where 



47r2 JJn (s — p^){s — p'^) 


(23) 


r d'^k 

pU = 3 ^y j^^{- 2 m)^ 6 {e) 6 [{k + p^k + p)^ - ml)] 
xTr[(/+ /: + m„) 7 a(i»+ /k + mb)ji,'y5{/i^ + md)jfj], 

Pia = 3i J -^^{-27ii)^6{k‘^)6[{k + p')^]S[{k + p)^-ml)] 

xTr[(/+ /: + m„) 7 „ 75 (j 6 + /c + mb)'ju75{A + (24) 
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and the corresponding integration region which can be solved via the Landan eqnation, 
is specihed by s' > 0 and s > ml + mis'/{ml — q^). Note that if {ml — q^)/q^ < s', there 
may be additional contribntions to the above integral becanse of pinching of the singnlarities 
on non-physical sheets |]^. However, at the hadron level, the contribntion from the higher 
states is approximated by the pertnrbative part, which starts from the thresholds Sq and 
Sq, we therefore stndy the behaviors of the form factors in the region: 0 < and 
{ml — g^)/g^ > Sg, where Sg is the threshold of the excited states, s' < Sg. In this region it is 
not necessary to worry abont possible contribntions from non-physical sheets. Considering 
only the leading power corrections of OPE of the correlation fnnction Eq. (14), we obtain 

13 




pert 

Pi 


{p,p';q^ 


+d2'^'gl{dd) 


Jo {s — p^){s'— p'^) 
2 / 2 


+ dQ^'gl{dd){bb) d 


+ + diidd) + diigsdaGd) 

P'gl{dd){uu) + ---, 


(25) 


where 



4^^572 ( 2 g'( 2 g - 2 ml - u)A - A[(mg - g^)^ - s'q'^] 
- md){2s' + 2ml - u)y 


(26) 


Ss' / 

f^aT = 4 ^ 2 y 7/2 (20AglTM(g ± s' - ml) - 2s ml] 

=FA[3Am — 2s'{{s — ml)"^ ± s'(s' -|- 2ml — g^))] 
-X\2s'±q^Tml)+C±), (27) 




3s' 


:(^X{3ml - 2q^ - s) -h 3A(3s' - s + g^)). 


(28) 


47r2A5/2 

with M = s-fs' —g^, A = —4ss', A = {s— ml){ml —q^) — s'ml, C+ = 4s'(20s'sA-|-6AA— 

As's), and C_ = 0. (See Ref. for the detailed calcnlation.) Note that in Eq. (25), the 
contribntion from the excited states has been approximately snbtracted dne to the npper 
integral limit, Sq or Sg, in the hrst term. To calcnlate the contribntion from the indnced 
condensate, we parametrize the bilocal condensate as in Ref. [|TH| : 















('l(0)S(0,x)(}(x)) = {qq)g(x\), 


( 29 ) 


with 


poo „ 

g{x\) = / daf{a)e-^E^/^, 

Jo 


(30) 


where the coordinates are Euclidean (x|; = —> 0) and g{x\) is the Euclidean space-time 
correlation function of the vacuum quarks. Our choice of the vacuum function is 

4 


/(«) = 


^-4a/m2 


rrtn 


(31) 


where ml{qq)= —{gsQf^^'^GfiuQ)- In this paper, we have adopted the convention = 9“ -|- 
igs^A^. If one adopts D" = 8°" — igs^A^, then ml{qq)= {gsqcr^'^G^uq)- The corresponding 
correlation function g{x‘^) of the vacuum quarks is of the monopole form 

1 


9{xe) = 


(32) 


1 -I- moxl/lG’ 

Since has been assumed to be spacelike, we have x\ — — 0- This choice leads to the 

empirical sea-quark distribution []^. By using the bilocal condensate parametrized above 
to Eq. (9), we obtain the relevant induced condensate in our study as 

(&(0)^(0,x)7^75m(x))^ = c?V*'^^(0|T(6(0)7^75M(x),h(|/)7^756(l/))|0) 

—the perturbative term) 


= AJ e*‘'"Xn,x2)dM, 


(33) 


with 


(j){u, x'^) = mb{uu) [ e u)^ 

Jo jj 


(34) 


In the following calculation, we adopt the approximation 1 since x^ ~ 0. Thus 

we obtain the contribution of the induced condensate, as shown in Fig. 2(b), 

U(j){u) 


n 


induced 


= du 


Aq,F^am- J, ““(pZ + X^’ 


(35) 
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and 2 {h) = n^'^“ 2 (b) = ^\^Pig%h) = 0) where (j){u) =4>{u,x‘^ = 0). After performing 

the donble Borel transform, Eq. (35) becomes 


B \T-rinduced 1 _ -rj 

[^^Ai,Fig.2{b)\ — -t* 


du 


u(l){u) 


u{p' + qY + (1 — u)p'‘^ + q^u{u — 1 )- 

X c- 


X 


M2 + M'2^ (M2 + M'2)2 


Note that the contribntion of the indnced condensate in Eq. (35) was ignored in Ref. [W 
The contribntion of indnced condensates as shown in Fig. 2(c) is given by 

■ninduced / / 2n ‘^dsi^d) A r 2\P^ - P'"^ - 

nx“2 (c)(p,p;?) = — 


9 


(p2 _ ’ 


Y^induced r /. 2\ _ 2g^Ydd) A( 2 ^ 

^^A2,Fig.2(c)lP)P ) ? ) — -^- XbuV^l J 

Yfinduced ( _ ‘^9s{dd) a ,2\ 

^^A,Fig.2ic)\PyP ^9 ) — ^ Xbu\9 J 


(p2 _ 77^2 ^p'4’ 


Tjinduced / /. 2\ _ 

^ ^V,Fig.2(c) [P J P ) — 


9 ^ (p 2 — 7772^p'4’ 

^9s{dd) y . 1 


9 


-xUq^) 


(p2 _ 


2'lr,'4’ 


(36) 


where we have approximated the non-local indnced condensate as a local one and the deh- 
nition of is 

(6(0)cr„^n(0))v 

= d^ye'^'^'‘^{ 0 \T{b{ 0 )ao,i 3 u{ 0 ),u{y)'jub{y))\ 0 ) — the pertnrbative term^ 

= -i'^''xLi9^)i9^aqy - dyuQa), 

(6(0)a„^75M(0))^ 

= d"^?/e*™( 0 |T( 6 ( 0 )(jQ,^ 75 n( 0 ), h(j/) 7 y 755 (|/))| 0 ) — the pertnrbative term^ 

= -i^’'xt{(f){9uaq0 - 9yuqa), (37) 

Similar to the case of two-point snm rnles, we can write down as a dispersion relation 

form, 

/*00 rj c 

X^^^\qY = / -^(Im 7r(s) - Im 7rP"’'*(s)), 

Jo s — q^ 
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and further adopt the simple model; 


Im 7r(s) = — rnjviA)) + Ini 7r^®^*(s)6'(s — 


The hnal values of nre 

fViA) 


V(A) _ 

Xbu _ 2 


r^o 


m 


fV(A) 


- q 


— — ds- 


^ Svr Jmi 


mb 
s — 


,ml 


m 


[l-2(—) + (—tl, 


(38) 


where /^=0.28 GeV^, /'^=0.67 GeV^, mjv=5.33 GeV, mj^=5.71 GeV, Sg=33 GeV^, and 
Sff=38.2 GeV^. The detailed calculation is shown in Ref. |^. Note that the uncertainties 
of our sum rule results are well controlled since in the following numerical analysis we hnd 
that the contribution of the induced condensate in Eq. (35) is about 10% to the Ai sum 
rule, and that the contribution of Eq. (36) to relevant sum rules is less than 4%. The rest of 
the results for di is collected in the Appendix. And dg, ds, d^\ d^\ and d® are represented 
pictorially by Figs. 3(d), 3(e), 3(f), 3(g), and 3(h), respectively. 

After equating the hadronic side to the quark side and then performing the double Borel 
transform B[/] on / in both variables and we derive the relevant form factor QGD 
sum rules as follows: 

fBml 

- -imB + 7 

9p 

Isml 1 


Qp ruB + mp 

1^2,n 

9p 

fsml 2 

Qp ruB + mp 




(39) 



(40) 


= B[nr'’1, 

(41) 



(42) 


III. RESULTS AND DISCUSSION 


A. The determination of Borel windows 


In numerical analyses of form factors, people usually assume, according to the empirical 
observation , that the suitable region of sum rules is at values of Borel masses twice as 
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large as that in the corresponding two-point functions |]l|,^ and thus extract results from 
larger Borel mass regions. However, we will provide a detailed analysis to show that such 
choice is not always correct. 


To obtain reliable Borel windows in which the form factor sum rules may be safely used 
in the analysis, we apply the differential operator —M'^d to both sides of Eq. (39) 

and then obtain the p mass sum rule, which is free of the parameters Pp, fs, and the form 
factors (Hi, A2, A, or V). This procedure is usually used in analyzing the mass sum rules in 
the two-point Green’s function approach [^,^. Analogously, we also obtain the p meson 
mass rules from Eqs. (40), (41), and (42). On the other hand, when applying the differential 
operator —M^d to both sides of Eqs. (39-42), we obtain four B meson mass sum 

rules. Totally, we obtain eight mass rules (four for the p meson and four for the B meson). 
In the numerical analysis, we use the the following set of parameters: 

rrib = 4.7 GeV, = rud = 0, (uu) = (dd) = —(240 MeV)^, = 0.6 — 0.8 GeV^, (43) 

where mb is a pole mass. The heavy quark condensate will be dismissed through this paper. 
To further improve the validity of the derived QGD sum rules, we have performed the 
following replacements: 

ml{qq) ml{qq)L-^, (qq) {qq)L^, (44) 


with 


L = 


ln((M'2M2)V2/A2 

ln(/i2/A2) 


(45) 


A=100 MeV, /i=500 MeV, and b = 11 —2n//3 with Uf being the number of unfrozen quark 
flavors. Note also that region R~^ << « ml, where R is the conhnement radius, exists 

an anomalous dimension 2/6 for both of the current operators b^u'^^q and b'y^q. In 

our numerical analysis of B meson decay form factors, the sum rules are studied at ~ 
VM'2M2 ^ 3.4 GeV2 [ see below, Eq.(47)]. Therefore, to obtain physical form factors, which 
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are independent of the subtraction scale, we put the factor /A^)] 

in the relevant form factor sum rules (the right hand side of Eqs. (39-42)). The requirement 
of obtaining reasonable mass results from these eight mass sum rules puts severe constraints 
on the choices of the parameters, Sq and Sg. With a careful study, we hnd that the best 
choice of parameters in our analysis is 

= 33 - 35 GeV^ sg = 1.3 - 1.4 GeVA (46) 


Gonsequently we hnd that 


1.5 GeV" < M" < 11.5 GeV", 0.9 GeV" < M'^ < 1.2 GeV' 


(47) 


are the reasonable choices through this paper. The resultant masses are mp=0.787±0.065 
GeV and mB=5.11±0.18 GeV. The detailed results are collected in Table I. As examples 
in Fig. 3(a) and Fig. 3(b), we plot the p mass and B mass (extracted from Fq.(39)), 
respectively, as a function of and Obviously, the study of both of the p and B mass 
sum rules is a gauge to understand the reliability of performing further numerical analyses 
on the form factors. In concluding this subsection, we will to discuss two “traditional” sum 
rule calculations by considering a three-point correlation function in the existing literature: 

The hrst is the work in Ref. [Q in which the authors use the vector current and the 
pseudoscalar current h'^^d as the interpolating helds for the p meson and for the B meson, 
respectively. For the form factor Ai, they obtain the following sum rule: 

^ mbPp 

3 f^o , Bo ^ j^_(s/M'^) -{s'(mbS -|-m,„(s — m^) 

V 2AV2 


, ds i ds e~ 
47r2 Jo Jo 


mbs'[{s — — ml) + s'mf] 

2/-1 , i2ml + 3mbmu - 2q^){ml - q^) 

- 

2ml + 3mbmu — 2q^ 3ml + 9mbmu — 4:q^ 


12M'2 


12M2 


, _^ml + 2mbmu-q‘^\^_mi/M^ 
+ ' 8M4 r 
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1 


1 


16as7r(gg)^mf, / 8 


9 


+ 


+ 


3ml 


2g2 


mt 


q 


V9M'2 9M2 iml-q^) 72M^ 


6M2M'2 


ml{ml — q^) {ml — q 


. 2\2 


+ 




( 48 ) 


72M6 ' 36M4M'2 

where A is defined as before. From their Ai sum rule, one can extract, following the procedure 
shown as above, the p meson mass sum rule. The result is illustrated in Fig. 4(a), from 
which we see that their result of p mass sum rule is not so stable as to determine the suitable 
Borel windows. Moreover, one can easily find that once a Borel window is determined, the 
resultant p mass from their A 2 or V sum rule is twice as large as that from the Ai sum rule. 
Note that if we adopt the varying external field to do the calculation again, we find that 
various induced condensates may enter complicatedly the sum rules for Ai, A 2 , and V and 
the contribution from induced condensates becomes a little big. Thus the estimate of the 
sum rule is less reliable. 

The second calculation is done by Ball and Braun [^. They use the tensor current 
and the pseudoscalar current b'y^d as the interpolating fields for the p meson and for 
the B meson, respectively. For the form factor ^ 2 , they have 


A, 


fBfpm^ 


B 


mb{mB + rup) 




Jo 


ds / ds e 


° ~ 'm\)u - 2ss' 


A3/2 


(s — mi){s — mg + 2s')vr — 3s'[(s — mly + 2(s — mys + ss'ju 
2ss'[(s — mlY + 2(s — ml)s' + 3ss'] 


A5/2 


mbml{qq) 2 jm2 
6M2M'2 
lQasTi{qq)‘^ / 1 


mt 


ml-q^ 


9 V6M'4 ^ 3M2M'2 ^ 36M4M'2 18M2M'4/ 


(49) 


where u and A are defined as before. Note that, for this A 2 sum rules, one obtains the same 
result if adopting the varying external field approach. In Fig. 4(b) we show the p mass 
sum rule extracted from Eq. (49). From Fig. 4(b) we obtain m^ < 0. The result seems to 
indicate that the tensor current da'^^u is not a good interpolating field for the p meson as 


mentioned in Ref. In the same reference, the authors have pointed out that the tensor 
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current cannot easily produce a stable p mass sum rule if the vacuum saturation hypothesis 
goes in the opposite direction: 


{qiiX‘qcni\'‘q) = -g(l +ft){?g)^ 


(50) 


with /ds —2. Moreover, the authors have also shown that the threshold of the excited 
states in their p meson mass sum rule turns out to be unphysically low, the sum rule is 
dominated by the continuum contribution, and the power corrections are rather large. 


B. Numerical analysis of form factors 


To study numerically the form factors, we adopt the set of parameters as shown in Eqs. 


(40), (43), and 


m 


/b = 160 MeV, Jp = 3.84. (51) 

The working Borel windows, which have been determined previously (Eq. (47)), are 8.5 
GeV^ < < 11.5 GeV^ and 0.9 GeV^ < < 1.2 GeV^. In these Borel ranges, the form 

factors are dominated by the leading perturbative bare loop; for the Ai form factor, the 
absolute value of the contribution of {ggdaGd) is less than 50% of the bare loop. Moreover, 
the contribution of induced condensate amounts to about 14% to the Ai sum rule, while 
less than 4% to 442, A, or V. Note that unlike light-cone sum rules §0. we cannot apply 
“directly” the heavy quark limit, rrih —>■ cx), to these sum rules since in that limit the series 
of the operator product expansion (OPE) may become unconvergent. But in the case of the 
B or D meson, our results indicate that the series of the OPE is in good convergence. In 
Fig. 3(a-d), we plot the Ai, A2, A, and V form factors at = 0, respectively, as a function 
of and M'^. We thus obtain the results on the form factors at = 0: 


44i(0)= 0.12 ±0.01, 

712 ( 0 )= 0.12 ±0.01, 

71(0)= 0.015 ±0.02, 

1/(0)= 0.15 ±0.02, (52) 
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where the error comes from the variation in the Borel mass, s^, Sq, or ttIq. The resultant 
values of form factors are half as large as the light-cone sum rule results or lattice 

QCD calculation We now consider the dependence of the various form factors. The 
variation of the form factors with is of great interest, since it probes the effects of strong 
interactions on the decay. As the property of discontinuity in Ref. [Q is mentioned (see also 
the discussion below Eq. (24)), the sum rules work well in the region (m^ — q^Y/q^ > Sq. 
Therefore, we could obtain the g^ dependence of the form factors over a wide range of g^ 
(from g^=0 up about to 9 GeV^). The g^ dependence of our form factors is given by 


F(g^) = F(0)(l-g7m^) 


2 \—n 


(63) 


where n = 1 for Ai, n = 2 for A 2 , A, and V, and the htted pole masses are m^^=5.45 
GeV, m^2=6.14 GeV, mA=5.98 GeV and my=5.78 GeV, respectively. Here the results 
are evaluated at the central values of the Borel mass ranges in Eq.(47). We hnd that 
our g^ dependence of the form factors is well consistent with the pole model ansatz by 
Kbrner and Schuler P] and recent lattice results [Q as well. In the following calculation, 
we will extrapolate our g^ dependence of form factors to all possible kinematic region. The 
pole model ansatz may be a good approximation for the form factor behavior since it is 
consistent with this sum rule calculation in the region: 0 GeV^ < g^ < 9 GeV^ and also in 
good accordance with the QGD power counting rules at large —g^ (the hard rescattering 
region). Moreover, by neglecting the light meson mass, we roughly obtain from Eq. (53) the 
relation F(g^)/F(0) ~ m^, where g^ = (m^ — rup)^. Therefore, our results agree with the 
prediction of heavy quark symmetry in the kinematic region near zero recoil (g^ ~ g^). 


a+ -|- a_~ m 


-3/2 

M ) 

-1/2 


~CIj^ Q— 

f 


m 


m 


-1/2 

M ) 

1/2 

M • 


(54) 


In Fig. 6(a) we plot the lepton-pair invariant mass spectrum dV/dq^ of the B ^ p i Ui 
decay together with dT^/dq"^, dT^/dq'^, and dV_/dq^. The solid curve is for dT/dq"^, and the 
long-dashed curve is for dTi/dq"^, the portion of the rate with a longitudinal polarized p in 
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the final state, the short-dashed curve is for dT^/dq^, the portion of the rate with a helicity 
minus p in the final state, while the dotted curve is for dVj^/dq^, the portion of the rate with 
a helicity positive p in the final state. Similarly, in Fig. 6(b) we plot the electron spectrum 
dr /dEe of the B ^ p i ui together with dTi/dEe-, dTj^/dE^^ and dT^/dEg. The solid curve 
is for dr /dEf., the long-dashed curve for dTi/dEg, the short-dashed curve for dT_/dEe, and 
the dotted curve for dT^/dEe- Both of the Figs. 6(a) and 6(b) are plotted in the B meson 
rest frame. 

From Fig. 6 we obtain dF/dg^ ~ dF i/dq^ and dF/ dE^, dF i/dEe at the maximum recoil 
region (g^ 0). This is quite reasonable because at low g^, the electron and antineutrino 

are nearly collinear, so that their net spin along their motion is zero. Since the B meson has 
spin zero, the energetic recoiling p meson must also have zero helicity. The helicity minus 
contribution is more weighted towards the large g^ value than the helicity zero contribution 
around the small g^ region. Our results also show that dF_/dg^ S> dF+/dg^ as expected 
from the left-chiral bi —> ul transition. Our results for the decay rate are given by 

F(S° ^ p^rPi)= (5.1 ± 1.0) X iKfeP X 10^2s-\ 

Tl/Tt= 0.85 ± 0.03, F+/F_ = 0.077 ± 0.012, 

T{B^ p+T-pr)= (3.1 ± 0.06) X iK^r X lO^V^ (55) 


Since the induced condensate like Eq. (35) does not contribute to the sum rules in Ref. 
for B^ — > therefore we take T{B^ —>■ 7r’''e“Pe)= (5.4 ± 1.6) x |FubP x 10^^s“^ 


and F(R° —vr+r i't)= (2.7 ± 0.07) x |146p x lO^^s ^ from Ref. |^. We obtain the ratios 
F(R° ^ p+e-pe)/ r(R0 ^ TT+e-pe)^ 0.94 and F(R0 ^ p+r-p^)/ T{B° vr+r-p^) ^ 1.15. 
Our result for the ratio F(R° —>■ p+e^Pg)/ F(R° —>■ vr+e^Pg) is a bit smaller than, but still 
consistent with the CLEO experimental value of 1.4±0.6. However, one should note 

that in the CLEO experiment, the reconstruction of the relevant events is model dependent. 
For instance, if the ISGW II model is used, the ratio F(.B° —>• p+e“Pe)/ T{B^ —>■ n~^e~Pe) is 
1.1T0.7. 

In closing this subsection, we consider the case of the D meson decays. By following the 
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same procedure as the case of —>• p+f' z/£, and using the parameters as in Ref. M 


we 


can obtain the same form factor results and their behaviors as in Ref. []^, except that 
(0) becomes 0.43± 0.04 and 74f^^*(0) becomes 0.61±0.04. The results for the form 
factors at = 0 are collected in Table II. The calculated decay rates read 


r(T)° ^ p~i^i^i)= 0.44 ± 0.08|14rfp X 

Tl/Tt= 0.58 ± 0.05, r+/r_ = 0.041 ± 0 . 002 , 

T{D+ K*^tui)= 0.51 ± O.OOIKsT X 

ri/rT= 0.99 ±0.06, r+/r_ = 0.19 ±0.02. (56) 


Taking 11481=0.975, we obtain B{D~^ =5.1±0.7 %. The experimental results 

H are B{D^ K*-i+ue) =4.8±0.4 %, Tl/Tt = 1.23 ± 0.13, r+/r_ = 0.16 ± 0.04. Our 
results are consistent with the existing experimental data, except that the value of r£,/rj’ is 
a little smaller than the experimental data. Further applications of this approach to various 


exclusive decay processes will be published elsewhere |31 


IV. SUMMARY 

In summary, we have used the varying external held approach of QCD sum rules to 
compute the form factors for the semileptonic decays B^ —> t~ vn. We have formulated 
this approach in a systematic way. By extracting both of the B meson and p meson mass 
sum rules, we can thus determine the reliable Borel windows in studying the relevant form 
factor sum rules. We also include induced condensate contributions, which have been ignored 
before, into the relevant sum rules. Therefore, we demonstrate that some QCD sum rule 
calculations in the literature are less reliable. Our results strongly support the pole model 
ansatz by Korner and Schuler on the dependence of the form factors. Combining with 
the previous analysis in Ref. |^, we obtain the ratio —>■ p+e“Pe)/ r(.B° —> 7r+e“Pe)~ 

0.94 and r(R° —p+r“pT-)/ r(i?° —>• ~ 1.15. 
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APPENDIX A: 


The quantities ds, ds, dg^^ df\ and d® dehned in Eq. (25) read 


4' 


2mu + TUd _md_ rndjml - q^) 

2p'2(p2 _ 2(p2 _ 77^2^2 2p'2(p2 _ 77^2^2 ’ 


(Al) 


dt 


1 r 6mdml 

12 I- (p2 — 


2md 

(p2 — 


+ 


4mfe — 2mu + Smd 


p'2(p2 



2mu + rrid 

p/4(p2 _ ^ 2 ^ 


2{md +rnu){jnl — q^) 4:{2ml — q'^)md + Gniuml 2md{ml — q'^Y 


p/4(p2 _ .^ 2^2 

(m^ — g^) ■ 
p'‘^{p^—m'lY -I’ 


pi2(jp. 


m 


2'|3 

b) 


pM(p2 


m 


2'|3 

b) 


(A2) 


d 


(1)^1 

6 


1 r 3 2 2 

~*”81 — ml) — ml)‘^ — ml)^ 

10ml + ‘28mbmu — 2q^ 6{ml — g^) -|- 7mbmu 0ml 

p'‘2(^p2 _ 7772^3 A p/4('p2 _ 7772^2 A ^^2 _ 7772)4 

2{ml — qYi^b + ^b^u — (f) Omliml + mbmu — g^)- 
A p"^{p‘^ — ml)^ p'‘^{p‘^ — ml)'^ -I’ 


(A3) 


d 


(2)Ai 

6 


2 r 1 1 

9 — m^) 


2muml + ml — 
p‘^{p^ — mDp'"^ 


Am,,ml 


p2(^p2 _ yjj2Ypl2 


(A4) 


pi3)Ai 

Ug 




4 1 


9p/4(p2 _ .^2) ’ 

rrid 

2p'2(p2 _ .^2)2’ 


(A5) 

(A6) 


19 





mu 


md 


6p"^{p‘^ — mlY Qpf‘^{jP‘ — miy 


(A7) 


j(l )^2 

ttg 


1 r 


+ 


27 Lj 9 ' 2 ('p 2 _ 77^2^3 p'4('p2 _ ,^^ 2 ^ 2 ] 


(A 8 ) 


d 


{2)A2 

6 


mr 


9ml Lp' 2 (^p 2 _ ^ 2 ^ 


pl2p2 


p/2 pi 2 


(A9) 


^(3)^2 

“6 


0 , 


(AlO) 




A 

3 


md 

2p'2(p2 _ jyj2y^ 


(All) 




A 

5 


6p"^{p‘^ — miy 


md 

2p/2(p2 _ 77^2^3 ’ 


(A 12 ) 


Qq 


1 r 2 3 1 

81 l-p' 2 (^p 2 _ ^2^3 p'4('p2 _ 77 ^ 2 ^ 2 ) ’ 


(A13) 




2 r 1 

9mf Lp'2('p2 _ ^ 2 ^ 



p/2p2 p'2p4 J ’ 


(A14) 


ci 


(3)A 

6 


0 , 


(A15) 


4 


md 

pl2(j)2 _ 77 ^ 2 ^ 2 ’ 


(A16) 


4 = 


4 )^ = - 

o 01 


m,, 




mdml 


md{ml - 


d‘i"'= * 


9ml 


2 p' 2 (p 2 — miy 

6p/4(p2 _ ^2y 

6ml 

2{ml - q^) 

p' 2 (p 2 _ 7772^4 

pli(j)2 — ffilY 

1 ml - 


1 



(A17) 


(A18) 


(A19) 


d 


(3)V 

6 


0 . 


(A 20 ) 


20 





REFERENCES 


[1] P. Ball, V. M. Braun, and H. G. Dosch, Pliys. Rev. D44, 3567 (1991); P. Ball, ibid, 48, 
3190 (1993). 

[2] V. M. Belyaev, A. Khodjamiryan, and R. Ruckl, Z. Phys. C60, 349 (1993). 

[3] A. Ali, V. M. Braun, and H. Simma, Z. Phys. C63, 437 (1993). 

[4] J. M. Flynn et al. (UKQCD coll.), Nucl. Phys. B461, 327 (1996); ibid, B476, 313 
(1996). 

[5] J.G. Kbrner and G. A. Schuler, Z. Phys. G46, 93 (1990). 

[6] M. Wirbel, B. Stech, and M. Bauer, Z. Phys. G29, 637 (1985). 

[7] N. Isgur, D. Scora, B. Grinstein, and M. Wise, Phys. Rev. D39, 799 (1989). 

[8] D. Scora and N. Isgur, Phys. Rev. D52, 2783 (1995). 

[9] H.-Y. Gheng, G.-Y. Gheung, and G.-W. Hwang, Phys. Rev. D55, 1559 (1997). 

[10] K.-G. Yang and W-Y. P. Hwang, Z. Phys. G73, 275 (1997). 

[11] P. Ball and V. M. Braun, Phys. Rev. D55, 5561 (1997). 

[12] H.-n. Li and H.-L. Yu, Phys. Rev. Lett. 74, 4388 (1995); H.-n. Li and H.-L. Yu, Phys. 
Lett. B353, 301 (1995). 

[13] A. V. Radyushkin and R. Ruskov, Nucl. Phys. B481, 625 (1996). 

[14] A. V. Smilga and M. A. Shifman, Yad. Fiz. 37, 1613 (1983) [ Sov. J. Nucl. Phys. 37, 
958 (1983)]. 

[15] K.-G. Yang, W-Y. P. Hwang, E. M. Henley, and L. S. Kisslinger, Phys. Rev. D47, 3001 
(1993). 

[16] R. E. Gutkosky, J. Math. Phys. 1, 429 (1960). 


21 



[17] One should note that not all tensor structures are independent. In the case of deriving 
the form factor g, we choose P^,T^a-PaT^,y, p'^T^a-p'aT^,u, PuTa/,, and 

as independent structures with = eaf.gpp'^p'', Sapu = ^.apupP^, and = e^pypp'^■ 
See also Ref. Hi- 


fis] H. Jung and L. S. Kisslinger, Nucl. Phys. A586, 682 (1995). The analytic form 
of the non-local condensate is uncertain. If a simple Gaussian-type mode for the 
coordinate dependence of the non-local condensate is adopted f{a)=5{a — tuq/A), 
(g( 0 )g(a;)) = (g( 0 )g( 0 )) xe“*^™'oA6 y_ Mikhailov and A. V. Radyushkin, Yad. Fiz. 
49,794 (1989) [Sov. J. Nucl. Phys. 49, 494 (1989)]), the contribution of this induced 
condensate to the Ai from factor is almost zero. However, Jung and Kisslinger have 
shown that if the simple Gaussian-type mode for the coordinate dependence of the non¬ 
local condensate is used, the resulting quark distribution is xq{x) oc 6{x — 8ml/M^), 
which is not acceptable(M^:the Borel mass). Moreover, since the induced condensates 
are not so important in my sum rules, the resultant uncertainties are thus reasonably 
small. 


[19] B. L. Ioffe and A. V. Smilga, Phys. Lett. B114, 353 (1982); B. L. Ioffe and A. V. Smilga, 
Nucl. Phys. B216, 373 (1983). 

[20] M. Neubert, Phys. Rev. D45, 2451 (1991). 

[21] M. A. Shifman, A. 1. Vainshtein, and V. 1. Zakharov, Nucl. Phys. B147, 385, 448 (1979); 
L. J. Reinders, H. Rubinstein, and S. Yazaki, Phys. Rep. 127, 1 (1985). 

[22] M. B. Voloshin and M. A. Shifman, Yad. Fiz. 45, 463 (1987) [Sov. J. Nucl. Phys. 45, 
292 (1987)]. 

[23] H. D. Politzer and M. B. Wise, Phys. Lett. B206, 681 (1988); ibid, 208, 504 (1988). 

[24] N. Isgur and M. B. Wise, Phys. Lett. B232, 113 (1989); 237, 527 (1990). 

[25] J. Govaerts, et ah, Nucl. Phys. B283, 706 (1987). 


22 



[26] S. J. Brodsky and G. P. Lepage, Phys. Rev. D22, 2157 (1980). 

[27] N. Isgur and M. B. Wise, Phys. Rev. D42, 2388 (1990). 


[28] M. Artnso, HEPSY 96-2, |hep-ex/96l00T]9 


[29] J. P. Alexander et al. (CLEO Collaboration), Phys. Rev. Lett. 77, 5000 (1996). 

[30] Particle Data Gronp, Phys. Rev. D54, 1 (1996). 


[31] K.-C. Yang, in preparation. 


23 



TABLES 


TABLE 1. Estimates of the p and B meson masses from Eqs. (39-42). 


mp (GeV) mB (GeV) 


Eq. (39) 

0.722~ 0.806 


5.15 ~ 5.28 

Eq. (40) 

0.820~ 0.844 


4.99 ~ 5.02 

Eq. (41) 

0.830~ 0.852 


4.92 ~ 4.95 

Eq. (42) 

0.764~ 0.803 


4.98 ~ 5.01 

Average 

0.787±0.064 


5.10±0.19 

TABLE II. Estimates of 

various weak semileptonic form factors for D - 

-P,K*. 

Ai(0) 

A2(0) 

A(0) 

V{0) 

^ p- 0.43±0.04 

0.57±0.08 

0.30±0.07 

0.98±0.11 

^ K*- 0.61±0.04 

0.67±0.08 

0.22±0.03 

l.lOiO.lO 
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Figure Captions 

Fig. 1. Diagrams for the propagator of Eq. (4). The heavy quark propagator is represented 
by the heavy line. 

Fig. 2. Diagrams for the correlation function of Eq. (20). The heavy quark propagator is 
represented by the heavy line. 

Fig. 3. The p mass and B mass, extracted from Eq. (39), plotted as a function of the square 
of the Borel masses and 


Fig. 4. The p mass as a function of the square of the Borel masses and M'^. In Fig. 
4(a) the p mass is extracted from Eq. (48), the Ai sum rule of |I|. In Fig. 4(b) the p mass 


is extracted from Eq. (49), the A 2 sum rule of jn 


Fig. 5. The ^i(O), ^ 2 ( 0 ), ^(0), and E(0) form factors plotted as a function of the square 
of the Borel masses and M'^. 


Fig. 6. (a) The lepton-pair invariant mass spectra dV/dq^ plotted as a function of q^. The 
solid curve is dT / dq^. The long-dashed curve stands for dF i/dq^, the portion of the rate with 
a longitudinal polarized p in the hnal state. The short-dashed curve stands for dV^/dq^, the 
portion of the rate with a helicity minus p in the final state, while the dotted curve is for 
(iF+/dg^, the portion of the rate with a helicity positive p in the hnal state, (b) The electron 
spectra dT/dEe of the B ^ p'^ i~ together with dT L/dE^, dT^/dE^, and dT_/dEe. The 
solid curve, the long-dashed curve, the short-dashed curve, and the dotted curve are for 
dT/dEe, dVi/dEe, dT_/dEe, and dT^/dEe-, respectively. 
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